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Faster Annealing Schedules for Quantum Annealing 
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New annealing schedules for quantum annealing are proposed based on the adiabatic the- 
] orem. These schedules exhibit faster decrease of the excitation probability than a linear 

schedule. To derive this conclusion, the asymptotic form of the excitation probability for 
| quantum annealing is explicitly obtained in the limit of long annealing time. Its first-order 

| term, which is inversely proportional to the square of the annealing time, is shown to be de- 

| termined only by the information at the initial and final times. Our annealing schedules make 

| it possible to drop this term, thus leading to a higher order (smaller) excitation probability. 

We verify these results by solving numerically the time-dependent Schrodinger equation for 
J> ' small size systems 

o ' 
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""^ ' 1. Introduction 

Oh, Quantum annealing (QA) is a novel generic algorithm to solve optimization problems. 1-3 

As can be seen from the fact that QA is also called quantum adiabatic evolution in the field of 
quantum computation, 4 QA is based on the adiabatic theorem of quantum mechanics. 5 One 



considers the time-dependent Hamiltonian consisting of two terms, the potential term and 
kinetic energy term, the latter being introduced by hand to solve the minimization problem of 
the former. The Hamiltonian initially has only the kinetic energy term and smoothly changes 
to eventually consist only of the potential term at the end of a process. If the Hamiltonian 
varies slowly, the adiabatic theorem guarantees that the system follows the ground state of the 
instantaneous Hamiltonian during the process. The idea of QA is that the non-trivial target, 
the ground state of the potential energy, is obtained after the adiabatic evolution starting 
from the trivial ground state of the kinetic energy. 

An important question is how the error rate of QA depends on the annealing time r. If 
r is too short, the change of the Hamiltonian is rapid and then the non-adiabatic transition 
may occur. Consequently the algorithm would miss the desired solution. An earlier discussion 
suggested that the residual energy defined by the energy difference between the solution 
obtained by an algorithm and true one is proportional to (logr) - ^ with £ ~ 6. 6 While this 
result is still under debate, 7 ' 8 the asymptotic behavior of the excitation probability in the large 
t limit is well-established. 8 According to the adiabatic theorem, 5 the excitation probability 
after the adiabatic evolution obeys the power law, r -2 . 
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To implement QA and derive the r -2 error rate, one ordinarily assumes that the Hamil- 
tonian depends linearly on time. However, there is no theoretical restriction on the annealing 
schedule, the changing rate of the Hamiltonian. In the present paper, we derive the explicit 
expression of the coefficient of the r -2 term in the excitation probability without this restric- 
tion of linearity. It is surprising that its upper bound is determined only by the information 
at the initial and final times. Based on this result, we find the condition that this coefficient 
vanishes and propose new faster annealing schedules having a much faster rate r _2m of error 
decrease. 

This article is organized as follows. In the next section, the upper bound of the excitation 
probability is derived for a general quantum system. We apply in §3 this result to QA and 
propose faster annealing schedules. In §4, we verify them by numerical simulations for the 
two-level system, the transverse field Ising model and the database search problem. The last 
section is devoted to conclusion. 

2. Upper bound for excitation probability 

Let us consider the general Hamiltonian which depends on time t only through the di- 
mensionless time s = t/r, 

Hit) = Hit/r) = H(s), (1) 
where r is the total annealing time. The state vector follows the Schrodinger equation, 

i±m)) = H(t)mt)h (2) 

or, in terms of the dimensionless time, 

i±\t{8))=TH{8)Ma)), (3) 

where we set h = 1. We assume that the initial state is chosen to be the ground state of H(0) 
and that the ground state of H(s) is not degenerate for < s < 1. Thus the state vector is 
expected to keep track of the instantaneous ground state of H(s) for sufficiently large r. The 
adiabatic theorem 5 provides a condition for the adiabatic evolution as 

r»max{A(s)}, (4) 



Aj( 8 ) 



(5) 



A; (s) 

where \j(s)) denotes the jth instantaneous eigenstate of H(s) with the eigenvalue £j(s) and 
Aj(s) = £j(s) — £o{s). We assume that |0(s)) is the ground state of H(s). 

Under the above adiabatic condition (4), the asymptotic expansion with respect to r yields 
the excitation amplitude at s = 1 as 

0(1)|^(1)) ~ ^ ds — (j(s)\ -Ai |0(*)> + 0(O. (6) 
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Using the integration by parts, we obtain the upper bound for the excitation probability as 



1 



0(i)IV(i)> < - 2 {Mo) + a 3 {i)Y + o(t-*). 



(7) 



This formula implies that the coefficient of the r -2 term is determined only by the state of 
the system at s = and 1 and vanishes if Aj(s) is zero at s = and 1. 

When the r~ 2 term vanishes, a similar calculation yields the next order term of the 
excitation probability. If H'(0) = H'(l) = 0, we obtain 



where we defined 



OXi)IV'(i)) 



< 



l (2), 



Af\0)+Af\l)} +0(t~ 



0(*)l— j-^w)> 



im+l 



ds T 



(8) 



(9) 



It is easy to see that the r -4 term also vanishes when H"(0) = H"(l) = 0. We generalize 
these results as follows: If the kth derivative of H(s) is equal to zero at s = and 1 for all 
k = 1, 2, ■ ■ ■ , m — 1, the excitation probability has the upper bound 



O'(i)hKi)) 

which is proportional to r~ 2m 



2 < — /4™)( )+4 m) i 



1 2 

ij +0(r- 



■2m-l> 



(10) 



3. Faster annealing schedules 

Although we have so far considered the general time-dependent Hamiltonian, the ordinary 
Hamiltonian for QA is composed of the potential term and the kinetic energy term, 

H(s) = {1 - f(s)} H khl + f(s)H pot . (11) 

The function f(s), representing the annealing schedule, satisfies /(0) = and /(l) = 1. Thus 
H(0) = i^kin and H(l) = H pot . The ground state of H pot corresponds to the optimal solution 
for the optimization problem. The kinetic energy is chosen so that its ground state is trivial. 
Therefore the above Hamiltonian connects the trivial initial state and the non-trivial desired 
solution. 

The condition for the r _2m term to exist is obtained straightforwardly from the results of 
the previous section because the Hamiltonian (11) depends on time only through the annealing 
schedule f(s). It is sufficient that the kth derivative of f(s) is zero at s = and 1 for k = 
1, 2, ■ ■ ■ ,m—l. We note that f(s) should belong to C m , that is, f(s) is an mth differentiable 
function whose mth derivative is continuous. 

Examples of the annealing schedules f m (s) with the r~ 2m error rate are the following 
polynomials: 

h(s)=s, (12) 



/ 2 ( s ) = s 2 (3-2s), 



(13) 
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/ 4 (s) = s 4 (35 - 84s + 70s 2 -20s 3 ). (15) 

The linear annealing schedule /i(s), which shows the r -2 error, has ordinarily been used in 
the past studies. 2,3 Although we here list only polynomials symmetrical with respect to the 
point s = 1/2, this is not essential. For example, f(s) = (1 — cos(7rs 2 ))/2 also has the r~ 4 
error rate because /'(0) = /'(l) = /"(0) = but f"(l) = -2vr 2 . 

4. Numerical results 

To confirm the upper bound for the excitation probability discussed in the previous section, 
it is instructive to study the two-level system with the Hamiltonian, 

H LZ (t) = -(l-f(tyj ha z - ao*, (16) 

where the a a (a = x, y, z) are the Pauli matrices. The energy gap of U\ji has the minimum 
2a at f(s) = 1/2. If the annealing time r is not large enough to satisfy eq. (4), the non- 
adiabatic transition occurs. The Landau-Zener theorem 9 ' 10 provides the excitation probability 

2 



P 

1 ex 



(i(i)IV'(i)) 



as 



ira 2 T 



where s* denotes a solution for f(s*) = 1/2. On the other hand, if r is sufficiently large, the 
system evolves adiabatically during the process. Thus the excitation probability has the upper 
bound (10), which is estimated as 

d m f 



Pcx ~ T 2m^2 +4a 2)m+2 



4h 2 a 2 r 



d m f 
J (0) 



ds 



+ 



-(1) 



(18) 



We numerically solved the Schrodinger equation (2) for this system (16) with the Flunge- 
Kutta method. Figure 1 shows the result for the excitation probability with annealing sched- 
ules (12)-(15). The initial state is the ground state of #lz(0). The parameters are chosen to 
be h = 2 and a = 0.2. The curved and straight lines show eqs. (17) and (18), respectively. In 
the small and large r regions, the excitation probability perfectly fits those two expressions. 

We carried out simulations of a rather large system, the Ising spin system with random 
interactions. The quantum fluctuations are introduced by the uniform transverse field. Thus, 
the potential and kinetic energy terms are defined by 

N 

tfpot = -£ J^to] -hJ2*!, (19) 

(ij) i=i 

N 

#km = -rj>f. (20) 
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r 



Fig. 1. The annealing-time dependence of the excitation probability for the two-level system (16) 
using schedules (12) to (15). The curved and straight lines show eqs. (17) and (18) for each 
annealing schedule, respectively. The parameters in eq. (16) are chosen to be h = 2 and a = 0.2. 



The initial state, namely the ground state of -ffkim is the all-up state along the x axis, 

\m) = ® lz+)i ^ z ~ h , (2i) 

i » 

where \z+)i and \z— )i stand for the eigenstates of af with eigenvalues 1 and —1, respectively. 

The residual energy E res , the energy difference between the solution obtained by the QA 
process and the exact one, is a useful measure for the error rate of QA. It has the same 
behavior as the excitation probability because it is rewritten as 

E ie8 = $(l)\H pot \j>(l)} - e (l) (22) 

= ^A J (l)|(j(l)|^(l))| 2 . (23) 

i>o 

Therefore E res is expected to be asymptotically in proportion to r _2m using faster annealing 
schedules. 

We investigated the two-dimensional square lattice of size 3x3. The quenched random 
coupling constants are chosen from the uniform distribution between —1 and +1, as 

shown in Fig. 2. The parameters are h = 0.1 and T = 1. Figure 3 shows the r dependence 
of the residual energy using annealing schedules (12)-(15). Straight lines representing r _2m 
(m = 1,2,3,4) are also shown for comparison. The data clearly indicates the r _2m -law for 
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-0.54 




Fig. 2. Configuration of random interactions {Jij} on the 3x3 square lattice which we investigated. 
The solid and dashed lines indicate ferromagnetic and antiferromagnetic interactions, respectively. 



large r. Note that the irregular behavior around E TCS ~ 10 comes from numerical rounding 
errors. 

Next, we apply the faster annealing schedule to the database search problem, finding an 
item in an unsorted database. Consider N items, among which one is marked. The goal of this 
problem is to find the marked item in a minimum time. The pioneering quantum algorithm 
proposed by Grover 11 solves this task in time of order y/~N, whereas the classical algorithm 
tests N/2 items on average. Farhi et al. 4 proposed a QA version of Grover 's algorithm and 
Roland and Cerf 12 found an annealing schedule which has the same complexity as Grover's 
algorithm. Although their schedule is optimal in the sense that the excitation probability by 
the adiabatic transition is equal to a small constant at each time, it has the r -2 error rate. 
We show that annealing schedules with the r~ 2m error rate can be constructed by a slight 
modification of their optimal schedule. 

Let us consider the Hilbert space which has the basis states \i) (i = 1, 2, ■ ■ ■ , N), and the 
marked state is denoted by \m). Suppose that we can construct the Hamiltonian (11) with 
two energy terms, 

H pot = l-\m){m\, (24) 
N N 

#kin = l-^£I»01- (25) 
i=l j=l 

The Hamiltonian H pot can be applied without the explicit knowledge of |m), which is the same 
assumption as in Grover's algorithm. The initial state is a superposition of all basis states, 

1 N 

^ N i=i 
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Fig. 3. The annealing-time dependence of the residual energy for the two-dimensional spin glass 
model with faster annealing schedules. The solid lines stand for functions proportional to T~ 2m 
(m = 1, 2, 3, 4). The parameters are h = 0.1 and T = 1. 



which does not depend on the marked state. The energy gap between the ground state and 
the first excited state, 



A 1 ( S ) = yi-4^/( S )(l-/( S )) ) (27) 

has a minimum at f(s) = 1/2. The highest eigenvalue £2(5) = 1 is (N — 2)-fold degenerate. 

To derive the optimal annealing schedule, we briefly review the results reported by Roland 
and Cerf. 12 When the energy gap is small, a non-adiabatic transition is likely to occur. Thus 
we need to change the Hamiltonian carefully. On the other hand, when the energy gap is not 
very small, too slow a change wastes time. Therefore the annealing schedule should be tuned 
to satisfy the adiabatic condition (4) in each infinitesimal time interval, that is, 

= S , (28) 

r 

where 5 is a small constant. In the database search problem, this condition is rewritten as 



^~ ld! S. (29) 



riVAi(s) 3 ds 

After integration under boundary conditions /(0) = and /(l) = 1, we obtain that 
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As plotted by a solid line in Fig. 4, this function changes most slowly when the energy gap 
takes a minimum value. It is noted that the annealing time is determined by the small constant 
5 as 



r=^^, (31) 



which means that the computation time is of order \/N similarly to Grover's algorithm. 

The optimal annealing schedule (30) shows the t~ 2 error rate because its derivative is 
not zero at s = and 1. It is easy to see from eq. (30) that a simple replacement of s with 
fm(s) fulfills the condition for the r~ 2m error rate. We carried out numerical simulations for 
N = 64 with such annealing schedules, f^is) = / op t (J m (s)), as plotted by dashed lines in 
Fig. 4. As shown in Fig. 5, the residual energy with f^(s) is proportional to r _2m . The 
characteristic time r c for the r~ 2m error rate to show up increases with m: Since the modified 
optimal schedule f^(s) has a steeper slope at s = 1/2 than f op t(s), a longer annealing time 
is necessary to satisfy the adiabatic condition (29). Nevertheless, the difference in slopes of 
fapt( s ) i s onr y a factor of 0(1), and therefore r c is still scaled as y/N. 

5. Conclusion 

In this paper, we have studied the error rate of QA based on the adiabatic theorem and 
proposed annealing schedules which present a faster decrease of the error rate with the anneal- 
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Fig. 5. The annealing-time dependence of the residual energy for the database search problem (N = 
64) with the optimal annealing schedules plotted in Fig. 4. The solid lines stand for functions 
proportional to r~ 2m (m — 1,2,3,4). 



ing time. The point is that the excitation probability is determined only by the information 
at t = and r when the system evolves adiabatically from the beginning to the end. Thus 
one can easily remove the usual r -2 error by tuning the annealing schedule. As examples of 
such schedules, we listed a few polynomials. 

Although our new faster schedules do not improve the characteristic time r c for the adia- 
batic evolution, these schedules can be used to produce a new faster schedule from the known 
efficient schedule with small r c . In general, it is difficult to obtain the optimal annealing 
schedule analytically. Nevertheless, if one once finds an efficient schedule for a specific prob- 
lem, schedules with smaller error rates are immediately constructed by a slight modification 
under the present scheme. 
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